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Abstract 


We  incorporated  the  CONWEP  blast  model  into  DYNA2D  and  DYNA3D.  It  works  as 
expected  and  appears  to  be  adequate  for  modeling  problems  such  as  vehicle  response  to  land 
mines.  The  model  accounts  for  the  angle  of  incidence  of  the  blast  wave,  but  it  does  not  account 
for  shadowing  by  intervening  objects  or  for  confinement  effects.  This  report  provides  FORTRAN 
hstings  and  directions  for  incorporating  the  model  in  DYNA2D,  DYNA3D,  and  associated 
preprocessors  and  postprocessors  and  suggests  changes  to  the  user  manuals. 
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1.  INTRODUCTION 


DYNA2D  [1]  and  DYNA3D  [2]  are  general  purpose  finite-element  programs  that  are  capable 
of  modeling  large  deformation  of  structures.  In  applications  such  as  blast  loading  of  the  floor  of  a 
vehicle  by  a  land  mine,  it  is  desirable  to  have  a  simple  blast-loading  model  rather  than  having  to 
explicitly  simulate  the  progress  of  the  shock  wave  fix)m  the  high  explosive  through  the  air  and  its 
interaction  with  the  structure. 

CONWEP  [3]  is  an  inplementation  of  the  empirical  blast  models  of  Kingery  and  Bulmash  [4, 5]. 
In  the  present  work,  the  CONWEP  model  was  incorporated  into  DYNA2D  and  DYNA3D  and  their 
interactive  graphics  preprocessors,  MAZE  [6]  and  INGRID  [7]. 

The  resulting  model  was  tested  by  simulating  a  blast-loading  experiment  reported  by  Pytleski  and 
Catherine  [8]  and  some  near-field  blast-impulse  experiments  by  Ewing  and  Huffington  [9].  The 
results  indicate  that  the  DYNA2D/CONWEP  and  DYNA3D/CONWEP  models  are  adequate  for  use 
in  engineering  studies  of  vehicle  response  to  the  blast  from  land  mines. 

2.  OBJECTIVES 

The  objective  of  this  work  was  to  provide  a  convenient  means  of  applying  blast  loading  on 
stmctures  in  the  DYNA2D  and  DYNA3D  codes,  without  having  to  ran  the  explosive  blast  problem 
exphcitly.  Analytic  formulas  fit  to  blast  data  were  available  from  the  work  of  Kingery  and 
Buhnash  [4,  5],  and  FORTRAN  subroutines  implementing  these  formulas  were  available  in  Hyde’s 
CONWEP  computer  program  [3]. 

3.  IMPLEMENTATION 

DYNA3D  already  has  a  blast-loading  capability  in  which  Erode  functions  can  be  applied  to 
structures  to  simulate  blast  from  a  nuclear  weapon.  We  extended  the  Erode  loading  function  to 
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incorporate  the  CONWEP  algorithms.  The  meaning  of  the  DYNA2D/DYNA3D  "IBRODE" 
parameter  was  extended  to  be: 

0:  no  blast  load 

1:  Erode  function 

2:  CONWEP  loading  function 

The  CONWEP  blast-toading  algorithms  were  extracted  from  the  CONWEP  package  and  modified 
where  necessary  to  conform  to  FORTRAN  77  syntax.  These  FORTRAN  subroutines  are  listed  in 
Appendix  A.  These  algorithms  provide  the  time  of  arrival,  peak  reflected  pressure,  peak  incident 
pressure,  and  duration  and  decay  of  the  incident  and  reflected  pressure. 

Input  data  required  by  the  CONWnP  model  are: 

•  weight:  equivalent  mass  of  TNT,  in  problem  mass  unit 

•  xO,  yO,  zO:  coordinates  of  the  point  of  explosion,  in  problem  length  units 

•  tO:  delay  time  between  when  the  DYNA  problem  starts  and  the  instant  of  explosion,  in  problem 
time  unit.  Can  be  negative. 

•  nunit:  units  switch: 

1:  pounds,  feet,  psi,  seconds 
2:  kilograms,  meters,  pascals,  seconds 
3:  dozens  of  slugs,  inches,  psi,  seconds 
4:  grams,  centimeters.  Megabars,  microseconds 
5:  conversion  factors  supphed  by  the  user 

•  isiuf:  surface  or  air  blast  switch 
1:  surface  blast 

2:  air  blast 
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In  addition,  the  DYNA2D/DYNA3D  model  requires  a  list  of  the  surface  segments  that  will 
expaience  the  blast  loading.  This  is  done  in  the  same  manner  as  with  the  Erode  model,  except  that 
the  load  curve  number  is  -2  instead  of  -1. 

The  Erode  model  does  not  account  for  the  angle  of  incidence  of  the  blast  wave  on  the  surface  of 
the  structure.  The  CONWEP  algorithms  do  account  for  angle  of  incidence  by  combining  the  reflected 
pressure  (normal-incidence)  value  and  the  incident  pressure  (side-on  incidence)  value.  Accordingly, 
we  modified  the  DYNA2D/DYNA3D  blast-loading  model  so  that  it  can  calculate  the  angle  of 
incidence  and  then  take  the  sum  • 

PressmeLoad  =  ReflectedPressure  •  cos^0 
IncidentPressure  •  (1  -i-  cos^0  -  2  cos0) 

When  COS0  is  negative  (i.e.,  the  surface  is  not  facing  the  point  of  explosion),  then 

PressureLoad  =  IncidentPressure, 

but  the  arrival  time  and  the  incident  pressure  are  not  adjusted  in  any  way  to  account  for  shadowing 
by  the  intervening  structure. 

The  DYNA2D/DYNA3D  blast-loading  model  does  not  take  into  account  any  confinement  or 
tunnel  effects  and  should  not  be  used  for  analyzing  such  problems. 

Except  for  the  addition  of  a  number  of  subroutines,  all  of  whose  names  begin  with  "conwep_", 
only  two  DYNA2D  subroutines  and  two  DYNA3D  subroutines  were  modified.  These  are 
DYNA2D's  subroutine  dynai  and  subroutine  fe2a,  and  DYNA3D's  subroutine  dynai  and  subroutine 
load.  These  modifications  are  shown  in  Appendices  E  and  C.  We  modified  DYNA2D's  graphics 
preprocessor,  MAZE,  to  read  and  write  CONWEP  data.  The  subroutines  that  were  modified  are 
blkdat,  meshol2,  and  meshoB,  and  we  added  suhaoutines  brodin  and  ndbrod.  The  changes  are  shown 
in  Appendix  D.  We  also  modified  DYNA3D's  graphics  preprocessor,  INGRID,  similarly.  The 
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modified  subroutines  are  dnopts,  tarodin,  and  ndbrod,  which  are  shown  in  Appendix  E.  Revisions  and 
additions  to  the  DYNA2D,  DYNA3D,  MAZE,  and  B^GRID  manuals  are  provided  in  Appendices  F, 
G,  H,  and  I,  respectively. 

4.  DISCUSSION 

4.1  Modeling  of  Blast-Loading  on  a  Heavy  Plate.  To  test  the  implementation  of  the  model  in 
DYNA3D,  we  ran  a  problem  in  which  a  2.32-kg  charge  was  detonated  as  a  surface  blast  at  a  standoff 
of  45.72  cm  firom  a  heavy  plate.  The  plate  was  modeled  as  a  very  dense  fluid  (20,000  g/cm^)  so  that 
its  acceleration  would  be  small  and  proportional  to  the  applied  pressure  load  and  that  its  velocity 
would  also  be  small  and  proportional  to  the  impulse.  Figure  1  shows  the  pressure  contours  on  the 
plate  at  1.1  ms  after  initiation  (multiply  the  velocity  by  1,000  to  determine  the  impulse  in 
Pascal-seconds).  Figure  2  shows  the  pressure-time  histories  at  10-cm  increments  of  radial  distance 
from  "ground  zero,"  and  Figure  3  shows  the  impulse-time  histories  at  the  same  locations.  The 
INGRID  input  file  for  this  problem  is  given  in  Appendix  J. 

4.2  Comparison  with  Experiments. 

4.2.1  Near-Field  Blast  Impulse  on  a  Heavy  Plate.  We  also  ran  a  series  of  problems  with 
DYNA2D,  using  experimental  conditions  described  in  Huffington  and  Ewing  [9].  They  measured 
the  refiected  inqjulse  near  spherical  charges  of  Pentolite  by  measuring  the  velocity  of  a  steel  plug  that 
was  accelerated  by  the  blast  impulse. 

In  DYNA2D,  the  table  was  modeled  as  a  very  dense  fluid  (10^  g/crr]?),  such  that  the  nodal 
acceleration  would  be  proportional  to  the  pressure  load  and  the  nodal  velocity  proportional  to  the 
impulse.  The  simulations  were  run  until  acceleration  was  complete  and  total  impulse  could  be 
determined.  A  MAZE  input  file  for  one  of  the  problems  is  given  in  Appendix  K. 
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The  initial  runs  revealed  a  problem  with  the  calculation  of  positive  phase  duration  when  the  scaled 
distance  was  less  than  the  range  of  applicability  claimed  by  Kingery  and  Bulmash  [4,  5].  Accordingly, 
we  modified  the  "conwep_tdur"  subroutine  to  use  a  constant  scaled  duration  in  such  cases. 

The  comparisons  between  the  DYNA2D/CONWEP  results  and  the  experimental  results  are 
shown  in  Table  1  and  Figure  4. 


c  Surface  blast  2.32kg  TNT  y=:45.72  cm 

TIME=  0.99999E-03 

CONTOURS  OF  MAX.  VELOCITY 

MIN=  0.142E+00  IN  ELEMENT  1 
MAX=  0.275E+01  IN  ELEMENT  615 


CONTOUR  VALUES 
A=:  O.OOE+00 
B=  2.50E-01 
C=  5.00E-01 
D=  7.50E-01 
E=1.00E+00 
F=1.25E+00 
G=  1.50E+00 
H=1.75E+00 
1=  2.00E+00 
J=  2.25E+00 
K=  2.50E+00 


Figure  1.  Contours  of  velocity  (proportional  to  impulse)  from  sample  problem. 
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Y  ACCELERATION 


c  Surface  blast  2.32kg  TNT  y=45.72  cm 


1.00E+04 


9.00E+03 


8.00E-I-03 


7.00E+03 


6.00E+03 


5.00E+03 


4.00E+03 


3.00E+03 


2.00E+03 


1.00E+03 


O.OOE+00 


O.OOE+00  1.50E-04  3.00E-04  4.50E-04  6.00E-04  7.50E-04  9.00E-04 


Sea 


Scaled  impulse  measured  by  Huffington  and  Ewing 


Figure  4.  Comparison  of  DYNA2D/CONWEP  model  with  experimental  results  bv  Huffington 
and  Ewing  TQI. 


Table  1.  Comparison  of  DYNA2D/CONWEP  Model  with  Experimental  Results 
by  Huffmgton  and  Ewing  [9]. 


Test 

Scaled  Distance 

ft/lb^^3 

Scaled  Impulse 
(DYNA2D/CONWEP) 

Scaled  Impulse 
[91 

T35 

0.48 

101. 

101. 

T66 

0.48 

100. 

98. 

T40 

0.38 

147. 

151. 

T72 

0.38 

149. 

150. 

T45 

0.29 

462. 

256. 

TlOO 

0.29 

149. 

152. 

T79 

0.29 

461. 

259. 

T107 

0.28 

467. 

274. 

T84 

0.19 

680. 

557. 

T95 

0.19 

679. 

588. 

T114 

0.10 

682. 

603. 

Till 

0.19 

683. 

579. 

T89 

0.16 

787. 

734. 

4.2.2  Dynamic  Deflection  of  an  Aluminum  Plate.  The  DYNA3D/CONWEP  model  was  "beta 
tested"  by  Ashutash  Jhaveri  and  John  Condon.  They  used  it  to  simulate  experimental  data  reported 
in  1993  by  Pytleski  and  Catherine  [8].  The  DYNA3D/CONWEP  predictions  were  in  excellent 
agreement  with  the  experimental  results.  Jhaveri  and  Condon's  results  were  presented  by  Condon, 
Gniazdowski,  and  Gregory  [10],  who  concluded  that  the  DYNA3D/CONWEP  model  is  sufficiently 
accurate  for  their  study  of  the  response  of  a  vehicle  floor  panel  to  blast  from  a  land  mine. 
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5.  CONCLUSIONS 


The  CONWEP  blast  model  has  been  incorporated  into  DYNA2D  and  DYNA3D.  It  works  as 
expected  and  appears  to  be  adequate  for  modeling  problems  such  as  vehicle  response  to  land  mines. 
The  model  accounts  for  the  angle  of  incidence  of  the  blast  wave,  but  it  does  not  account  for 
shadowing  by  intervening  objects  or  for  confinement  effects. 
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APPENDIX  A: 

CONWEP  SUBROUTINES,  CONVERTED  TO  FORTRAN  77 
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subroutine  conwep_angle(x  1  ,x2,x3,x4,y  1  ,y2,y3,y4, 

&  zl,z2,z3,z4,s,x5,y5,z5,cosa) 

c  Find  the  incidence  angle  from  xO,yO,zO  to  the  face  that 
c  has  comers  (xl,yl,zl) ...  (x4,y4,z4). 
c  Glenn  Randers-Pehrson,  22  April  1994 

c  "s"  is  the  load  curve  multiplication  factor.  If  negative, 
c  assume  that  the  nodes  are  numbered  in  reverse  order 

common  /conwep_input/  wtnt,xO,yO,zO,tO,isurf 

c  Find  vectors  connecting  opposite  pair  of  comers 
c  Call  the  vectors  "a"  and  "b" 

ax  =  x4  -  x2 
ay  =  y4  -  y2 
az  =  z4  -  z2 

bx  =  x3  -  xl 
by  =  y3  -  yl 
bz  =  z3  -  zl 

c  Take  cross  product 
c  Call  the  result  a  vector  "e" 

ex  =  by*az  -  ay*bz 
ey  =  bz*ax  -  az*bx 
ez  =  bx*ay  -  ax*by 

c  Normalize.  The  unit  vector  points  out  of  the  surface 

emag  =  SQRT(ex*ex  +  ey*ey  +  ez*ez) 
if(emag.eq.O.)emag=l . 
ex  =  ex/emag 
ey  =  ey/emag 
ez  =  ez/emag 

c  Find  centroid 

if(x3.eq.x4  .and.  y3.eq.y4  .and.  z3.eq.z4)then 
x5=(xl+x2+x3)/3. 
y5=(yl+y2+y3)/3. 
z5=(zl+z2+z3)/3. 
else 
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x5=(x  1  +x2+x3+x4)/4. 
y5=(y  1  +y2+y3+y4)/4. 
z5=(zl+z2+z3+z4)/4. 
endif 

c  Find  vector  from  centroid  of  face  to  center  of  blast 
c  Call  the  result  a  vector  "f ' 

fx  =  xO  -  x5 
fy  =  yO  -  y5 
fz  =  zO  -  z5 

c  Normalize.  The  unit  vector  points  toward  the  center  of  blast 

fmag  =  SQRT(fx*fx  +  fy*fy  +  fz*fz) 
if(fmag.eq.O.)fmag=l . 
fx  =  fx/fmag 
fy  =  fy/fmag 
fz  =  fz/fmag 

c  find  angle  between  vectors  "e"  and  "f 

edotf  =  ex*fx  +  ey*fy  +  ez*fz 
cosa  =  MIN(  MAX  (edotf, -1.0),  1.0  ) 
if  (s.le.O.)  cosa  =  -cosa 

end 


subroutine  conwep_blast(xl ,x2,x3,x4,yl ,y2,y3,y4, 

&  zl,z2,z3,z4,s,t,p) 

c 

common  /conwep_input/wtnt,xO,yO,zO,tO,isurf 
c 

common  /conwep_units/  iunit,uft,ulb,upsi,ums 
real* 8  a,  b 

logical  warn 
save  warn 
data  wam/.true./ 
c 

c  This  subroutine  invokes  Conwep  logic  to  compute  the  pressure 
c  load  P  on  a  surface  at  a  field  point  (x,y,z),  due  to  a  conventional 
c  TNT  charge  undergoing  a  hemispherical  surface  burst  (isurf=l)  or 
c  spherical  air  burst  (isurf=2)  located  at  point  (x0,y0,z0)  and 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C' 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


time  tO  relative  to  the  start  of  a  DYNA3D  simulation.  P  is 
returned  in  psi,  Pascals,  or  Megabars,  depending  on  the  value 
of  nunit. 

The  user  must  supply  the  following  as  input  to  DYNA3D,  subroutine 
conwep_blasti: 

weight  -  Equivalent  weight  or  mass  of  TNT,  in  lbs.  Kilograms, 
dozens  of  slugs,  or  grams,  respectively. 

(xO,yO,zO)  -  Coordinates  of  point  of  explosion  w.r.t.  global 
origin,  in  feet,  meters,  inches,  or  cm 
to  -  Delay  time,  if  any,  between  when  problem 

starts  and  instant  of  explosion,  seconds  or  microseconds 
nunit  -  Units  switch 

=  1 ,  engineering  units  (feet,  pounds,  psi,  seconds) 

=  2,  S.I.  units  (meters,  kilograms,  pascals,  seconds) 

=  3,  English  units  (inches,  dozens  of  slugs,  psi,  seconds) 

=  4,  metric  units  (cm,  grams.  Megabars,  microseconds) 

=  5,  conversion  factors  are  supplied  by  the  user 

Also,  to  activate  this  logic,  the  user  must  use  NL  =  -2  when 
defining  the  pressure  load  input  to  DYNA3D. 

The  equations  from  BRL  Technical  Report  ARBRL-TR-02555  are  used 
to  find  the  incident  pressure  P,  time  of  arrival  ta,  time  of  pulse 
duration  to,  impulse  impo,  peak  incident  overpressure  pso, 
reflected  impulse  impr,  and  peak  reflected  overpressure  pro  for 
a  field  point  at  slant  range  R  from  the  point  of  explosion 
(xO,yO,zO)  for  a  conventional  charge  of  TNT.  The  peak  incident 
(side-on)  and  peak  reflected  (head-on)  pressures,  together  with 
the  computed  incidence  angle  alpha,  are  used  to  determine  the 
pressure  load  on  the  face. 

Description  of  variables  (units  given  are  as  used  internally): 

a  -  decay  coefficient  of  incident  pressure  vs.  time  history, 
where  p(t)  =  pso*(l-t/to)*exp(-a*t/to) 

This  parameter  is  computed  in  function  decay, 
b  -  decay  coefficient  of  reflected  pressure  vs.  time  history, 
where  p(t)  =  pro*(l-t/to)*exp(-b*t/to) 

This  parameter  is  computed  in  function  decay, 
cosa  -  cos(incidence  angle)  (1  =  head  on,  0  =  side  on) 
isurf  -  1  =  surface  burst,  2  =  air  burst  (set  to  1  currently) 

P  -  incident  pressure  at  field  point  (x,y,z,t)  and 
time  t;  computed  in  routine  pressure,  psi 
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c  pro  -  peak  reflected  overpressure,  psi 
c  computed  in  routine  params  by  calling  function  pref 

c  pso  -  peak  incident  overpressure,  psi 
c  computed  in  routine  params  by  calling  function  pine 

c  R  -  slant  range,  =  sqrt((x-x0)**2+(y-y)**2+(z-z0)**2),  ft 
c  s  -  multiplication  factor.  If  negative,  assume  nodes  are 
c  numbered  in  reverse  order 

c  t  -  current  DYNA3D  time 
c  ta  '  arrival  time,  msec 

c  Computed  by  routine  params  by  use  of  function  tarr. 

c  to  -  positive  phase  duration  time,  msec 
c  Computed  by  routine  params  by  use  of  function  tdur. 

c  ts  -  shifted  time,  =  t  -  tO 
c  wtnt  -  equivalent  weight  of  TNT,  lb 
c  x,y,z  -  coordinates  of  field  point  w.r.t  global  origin  where 
c  incident  pressure  is  to  be  computed,  default  unit 

c  is  meters 

c  zlant  -  scaled  slant  range,  ft/lb**  1/3 
c  Computed  in  routine  params. 

c  zlog  -  logarithm  of  scaled  slant  range 
c  Computed  in  routine  params. 

c 

c  Compute  centroid  of  face  and  incidence  angle 
c 

call  conwep_angle(x  1  ,x2,x3 ,x4,y  1  ,y2,y3,y4,zl  ,z2,z3,z4, 

&  s,x,y,z,cosa) 

c 

c  Compute  slant  range  (distance  between  field  point  and  the  point  of 
c  explosion).  Convert  to  appropriate  units, 
c 

R  =  sqrt((x-x0)**2  +  (y-y0)**2  +  (z-z0)**2)*uft 
c 

c  Compute  pso,  ta,  to,  zlant,  and  zlog  for  this  point, 
c 

call  conwep_params(isurf,R,ta,to,pso,pro,wtnt,zlant,zlog,a,b) 
c 

c  Get  shifted  time  ts  and  express  it  in  terms  of  milliseconds  for 
c  internal  use. 
c 

ts  =  (t  -  tO)*ums 

if(  (ta  .It.  -tO*ums)  .and.  warn)  then 
write(*,*)'  Caution:  arrival  time  <  0  in  conwep_blast' 
write(*,*)'  ta  ■=',ta 
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if  (t.ne.O.)wam=.false. 
endif 
c 

c  Compute  the  incident  pressure  P(x,y,z,ts). 
c 

call  conwep_press(a,b,ts,P,ta,to,cosa,pro,pso) 
c 

c  Convert  pressure  to  appropriate  units  before  returning  to 
c  calling  program, 
c 

p  =  P/upsi 

return 

end 


SUBROUTINE  conwep_blasti 
C 

C  READ  AND  PRINT  BLAST  FUNCTION  DATA  FOR  USE  IN  SUBROUTINE 
C  CONWEP_BLAST  (CALLED  BY  SUBROUTINE  FE2A) 

c  weight  -  Equivalent  weight  or  mass  of  TNT,  in  lbs.  Kilograms, 
c  dozens  of  slugs,  or  grams,  respectively, 

c  (xO,yO,zO)  -  Coordinates  of  point  of  explosion  w.r.t.  global 
c  origin,  in  feet,  meters,  inches,  or  cm 

c  tO  -  Delay  time,  if  any,  between  when  problem 
c  starts  and  instant  of  explosion,  seconds  or  microseconds 

c  nunit  -  Units  switch  (default  2) 

c  =  1 ,  engineering  units  (feet,  pounds,  psi,  seconds) 

c  =2,  S.I.  units  (meters,  kilograms,  pascals,  seconds) 

c  =3,  English  units  (inches,  dozens  of  slugs,  psi,  seconds) 

c  =4,  metric  units  (cm,  grams.  Megabars,  microseconds) 

c  =5,  conversion  factors  are  supplied  by  the  user 

c  isurf  -  Surface  or  Air  Blast  switch  (default  2) 
c  =  1 ,  surface  blast 

c  =2,  air  blast 

c 

common  /conwep_input/wtnt,xO,yO,zO,tO,isurf 

common  /conwep_units/  iunit,uft,ulb,upsi,ums 

CHARACTER*80  TXTS,MSSG 
C 

CALL  GTTXSG  (TXTS,LCOUNT) 
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non 


READ  (UNIT=TXTS,FMT=100,ERR=90) 
1  WEIGHT, X0,Y0,Z0, TO, nunit,nsurf 


iunit  =  nunit 
if(iunit.eq.0)iunit=2 


isurf  =  nsurf 
if(isurf.eq.0)isurf=2 


if(iunit.eq.5)then 

CALL  GTTXSG  (TXTS,LCOUNT) 

READ  (UNIT=TXTS,FMT=  1 00,ERR=90)uft,ulb,upsi,ums 
else 

call  conwep_stunit 
endif 

c  convert  weight  to  appropriate  units  for  internal  use. 

wtnt  =  weight*ulb 
c 

WRITE(12,1 10)  WEIGHT, WTNT,X0,Y0,Z0,T0,IUNIT,ISURF 
WRITE(12,120)  ulb,uft,ums,upsi 
C 

RETURN 


TERMINATE  DUE  TO  BADLY  FORMATTED  DATA 


90  MSSG='  BLASTI:  ERROR  READING  BLAST  INPUT  DATA' 
CALL  TERMIN  (TXTS,MSSG,LCOUNT,l) 


100  FORMAT(5E10.0,2i5) 

1 10  FORMAT(//'  CONWEP  OVERPRESSURE', 

1  '  CALCULATION'// 

1  4X,'WEIGHT  ....',1PE1L4,'....  WEIGHT  IN  LB . ',1PE1L4// 

3  4X,'XB0 . M PEI  1.4// 

4  4X,'YB0 . ',1  PEI  1.4// 

5  4X,'ZB0 . ',1  PEI  1.4// 

6  4X,'TB0 . ',1PE11.4// 

7  4X,'IUNIT . ',11 1// 

7  4X,'ISURF . ',11 1//) 

1 20  fomiat( 


1  4X,'LB/MASS  UNIT . ',1PE11.4// 

2  4X,'FT/LENGTH  UNIT . ',1PE1 1.4// 

3  4X,'MILLISEC/TIME  UNIT . ',1PE1 1.4// 

4  4X,'PSFPRESSURE  UNIT . ',1PE1 1.4//) 

END 
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blockdata  conwep_brl 

common  /conwep_input/wmt,xO,yO,zO,tO,isurf 
common  /conwep_brlpso/  cpso(0: 1 1 ,2) 
common  /conwep_ambient/  cO,  pO,  gO,  rhoO 
common  /conwep_units/  iunit,uft,ulb,upsi,ums 
c 

c  Coefficients  for  incident  pressure  equations 
c 

data((cpso(i,j),i=0,ll),j=l,2)  / 

&  1.9422502013,  -1.6958988741, 

& -0.154159376846,  0.514060730593, 

&  0.0988534365274,  -0.293912623038, 

&  -0.0268112345019,  0.109097496421, 

&  0.001 628467563 1 1  ,-0.02 1 4631030242, 

&  0.0001456723382,  0.00167847752266, 

&  1.77284970457,  -1.69012801396, 

&  0.0080497359 1951,  0.336743 1 1 4941 , 

&  -0.005 1 622635 1 334,-0.08092286 19888, 

&.  -0.00478507266747,  0.00793030472242, 

&  0.0007684469735,  3*0.0/ 

data  cO/ 1 1 1 6./p0/ 1 4.696/gO/ 1 .4/rho0/0.0765/ 
end 


double  precision  function  conwep_decay(pO,iO,td) 
real*4  pO,iO,td 
real*  8  a,fa,fpa 
c 

c  find  rate  of  decay  for  pressure  assuming: 
c 

c  p(t)  =  pO  *  [  1  -  (t-ta)/td]  *  exp[-a*(t-ta)/td] 
c  (ffiedlander's  equation) 
c 

c  where  a  is  a  decay  coefficient,  integrating  this  equation 
c  over  time  gives  the  impulse: 
c 

c  iO  =  pO  *  td  *  [  a  +  exp(-a)  -  1  ]  /  a**2 
c 

c  find  f(a)  =  a**2  -  pO*td/iO*[  a  -i-  exp(-a)  -  1]  =  0 
c  or: 

c  a**2  /  [  a  -f  exp(-a)  -  1]  -  pO*td/iO  =  0 
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c 

c  for  large  a,  exp(-a)  approaches  0,  and 
c  a**2  /  (a  -  1)  approaches  a  +  1 
c 

c  initial  guess  a  =  p0*td/i0  -  1 
c 

ptoi  =  p0*td/i0 
c  initial  guess: 
a  =  ptoi  -  1 . 

1  fa  =  a*a  -  ptoi*(a  +  exp(-a)  -  1.) 
fpa  =  2*a  -  ptoi*(l.  -  exp(-a)) 
a  =  a  -  fa/fpa 

if(abs(fa)  .gt.  1  .e-6)  go  to  1 
conwep_decay  =  a 
return 
end 


subroutine  conwep_params(ISURP,R,ta,to,pso,pro,WTNT,z,zlog,a,b) 
c 

c  purpose:  use  the  equations  from  brl  technical  report 
c  arbrl-tr-02555  to  find  the  incident  pressure  and  time  of 
c  arrival  for  a  point  in  space  at  a  given  slant  range  from 
c  conventional  explosion, 
c 

c  description  of  variables: 
c  a  -  decay  coefficient  of  incident  pressure 
c  vs.  time  history,  where 

c  p(t)  =  pso*(l-t/to)*exp(-a*t/to) 

c  b  -  decay  coefficient  of  reflected  pressure 
c  vs.  time  history,  where 

c  p(t)  =  pro*(l-t/to)*exp(-b*t/to) 

c  pso  -  peak  incident  overpressure,  psi 

c  pro  -  peak  reflected  overpressure,  psi 

c  R  -  slant  range,  ft 
c  ISURF  -  1  =  surface  burst,  2  =  air  burst 
c  ta  -  arrival  time,  msec 
c  to  -  positive  phase  duration,  msec 
c  WTNT  -  equivalent  weight  of  tnt,  lb 
c  z  -  scaled  slant  range,  ft/lb**  1/3 
c  zlog  -  logarithm  of  scaled  slant  range 
c 


real* 8  a,  b,  conwep_decay 
real  impo,  impr 
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logical  warn 
save  warn 
data  wam/.true./ 

w3  =  WTNT**(l./3.) 
z  =  R  /  w3 
zlog  =  aloglO(z) 

if(ISURF.eq.l)  then 
zlo  =  0.45 
else 

zlo  =  0.37 
endif 

ta  =  conwep_tarr(ISURF,ZLOG)  *  w3 
to  =  conwep_tdur(ISURF,ZLOG)  *  w3 
impo  =  conwep_xinips(ISURF,ZLOG)  *  w3 
impr  =  conwep_xinipr(isurf,zlog)  *  w3 
pso  =  conwep_pinc(ISURF,ZLOG) 
pro  =  conwep_pref(isurf,zlog) 


if(z  .ge.  zlo)  then 

if(pso*to/inipo  .le.  2.5  .or.  pro*to/impr  .le.2.5)then 
if(wam)then 

write(*,*)  'conwep_paranis:  resetting  to  from  ',to 
to  =  2.5*max(impo/pso,impr/pro) 
write(*,*) '  to  ',to 

wam=.false. 
else 

to  =  2.5*max(impo/pso,inipr/pro) 
endif 
endif 

a  =  conwep_decay(pso,impo,to) 
b  =  conwep_decay(pro,impr,to) 
else 

if(pso*to/impo  .le.  2.  .or.  pro*to/impr  .le.2.)then 
if(wam)then 

write(*,*)  'conwep_params:  resetting  to  from  ',to 
to  =  2.0*max(impo/pso,inipr/pro) . 
write(*,*) '  to  ',to 

wam=.false. 
else 

to  =  2.0*max(impo/pso,inipr/pro) 
endif 
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endif 
a  =  0. 
b  =  0. 
endif 

return 

end 


real  function  conwep_pinc(isuif,zlog) 
c 

c  purpose;  find  the  incident  pressure  due  to  the  detonation 
c  of  a  1  lb  equivalent  tnt  charge,  equations  are  from  BRL 
c  technical  report  ARBRL-TR-02555.  pressure  is  returned  in 
c  units  of  psi. 
c 

c  description  of  variables; 
c  isurf  - 1  for  surface  burst,  2  for  air  burst 
c  zlog  -  logarithm  (base  10)  of  scaled  range 
c 

common  /conwep_brlpso/  cpso(0;  1 1 ,2) 
c 

u  =  -0.756579301809  +  1. 35034249993 *zlog 
conwep_pinc  =  cpso(l  1, isurf) 
do  10  i  =  10, 0,-1 

10  conwep_pinc  =  conwep_pinc*u  +  cpso(i,isurf) 
conwep_pinc  =  10.**conwep_pinc 
return 
end 


real  function  conwep_pref(isurf,zlog) 
c 

c  purpose;  find  the  normally  reflected  pressure  due  to  the  detonation 
c  of  a  1  lb  equivalent  mt  charge,  equations  are  from  BRL 
c  technical  report  ARBRL-TR-02555.  pressure  is  returned  in 
c  units  of  psi. 

c 

c  description  of  variables; 
c  isurf  - 1  for  surface  burst,  2  for  air  burst 
c  zlog  -  logarithm  (base  10)  of  scaled  range 
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parameter  (ns=l  1,  nf=9) 

real  csurf(0:ns),  cfree(0:nf) 

data  csurf  /  2.56431321138,  -2.21030870597, 

&  -0.218536586295,  0.895319589372, 

&  0.24989009775,  -0.569249436807, 

&  -0.11791682383,  0.224131161411, 

&  0.0245620259375,  -0.0455116002694, 

&  -0.00190930738887, 0.00361471 193389  / 

data  cfree  /  2.39106134946,  -2.21400538997, 

&  0.035119031446,  0.657599992109, 

&  0.0 1 41 8 1 895 1 887,  -0.24307663623 1 , 

&  -0.0158699803158,  0.0492741184234, 

&  0.00227639644004,-0.00397126276058  / 

if  (isurf.eq.l)  then 

u  =  -.789312405513  +  1. 3663771 9229*zlog 
pref  =  csurf(ns) 
do  10  i  =  ns-1,0,-1 
1 0  pref  =  pref*u  -t-  csurf(i) 
else 

u  =  -0.756579301809  +  1 .35034249993*zlog 
pref  =  cfree(nf) 
do  20i  =  nf-l,0,-l 
20  pref  =  pref*u  -i-  cfree(i) 
endif 

conwep_pref  =  10.**pref 

return 

end 


subroutine  conwep_press(a,b,ts,p,ta,to,cosa,pro,pso) 
real*  8  a,b 

common  /conwep_units/  iunit,uft,ulb,upsi,ums 
c 

c  Compute  pressure  at  shifted  time  ts  for  this  field  point  in  space, 
c 

c  The  shifted  time,  ts  (in  milliseconds),  is  defined  by 
c 

c  ts  =  (t  -  tO)*ums 

c 

c  where  t  =  current  DYNA3D  elapsed  time,  seconds 
c  to  =  delay  time  for  the  explosion,  seconds 
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c  to  >  0.0  means  explosion  takes  place  tO  seconds  after 
c  DYNA3D  calculation  starts 

c  to  =  0.0  means  explosion  coincides  with  start  of 

c  DYNA3D  calculation 

c  to  <  0.0  means  explosion  occurs  tO  seconds  before 
c  DYNA3D  calculation  starts, 

c 
c 

c  If  the  pressure  front  has  not  arrived  here  yet  ( ts  .LT.  ta ), 
c  and  this  will  include  the  case  when  ts  <  0.0  (i.e.  (t  -  tO)  <  0.0 ), 
c  then  the  pressure  at  this  field  point  is  zero, 
c 

c  Parameters  a,  ta,  to,  and  pmax  are  functions  of  the 
c  charge  parameters  and  the  slant  range  distance  from  this  field 
c  point  to  the  point  of  explosion.  These  parameters  must  have  been 
c  correctly  computed  before  entering  this  routine, 
c 

iffts.ge.ta)  then 

exa  =  exp(-a*(ts-ta)/to) 
exb  =  exp(-b*(ts-ta)/to) 
poscosa=MAX(cosa,0.0) 

p  =(pso*exa*(l.+poscosa-2.*poscosa**2)+pro*exb*poscosa**2) 
&  *(l-(ts-ta)/to) 

p  =  max(p,-14.7) 

else 

p  =0.0 

endif 

c 

return 

end 


SUBROUTINE  conwep_stunit 
common  /conwep_units/  iunit,uft,ulb,upsi,ums 
IF(IUNIT  .EQ.  0.  .or.  iunit  .eq.  1)  THEN 
c  Engineering:  ft-lb-sec-psi 
UFT  =  1. 

ULB  =  1. 

UPSI=  1. 
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urns  =  1000. 

ELSEIF  (lUNTT  .eq.  2)  THEN 
c  SI:  m-kg-sec-Pascal 
UFT  =  1.0/0.3048 
ULB  =  1.0/0.45359 
UPSI  =  1./6894.80 
urns  =  1000. 

ELSEIF  (lUNTT  .eq.  3)  THEN 
c  English:  inch-12slugs-seconds-psi 

UFT  =  1./12. 

ULB  =386.09 
UPSI  =  1. 
urns  =  1000. 

ELSEIF  (lUNIT  .eq.  4)  THEN 
c  metric:  cm-g-microsec-Mbar 

UFT  =  1.0/30.48 
ULB  =  .001/0.45359 
UPSI=  l./(6894.8*l.e-ll) 
urns  =  .001 
ENDIF 
RETURN 
END 


real  function  conwep_taiT(ISURF,ZLOG) 
c 

c  purpose:  find  the  scaled  time  of  arrival  for  the  detonation  of  a 
c  1  lb  equivalent  tnt  charge,  equations  are  from  BRL 
c  technical  report  ARBRL-TR-02555.  arrival  time  is  returned 

c  in  msec/lb**(l/3). 

c 

c  description  of  variables: 
c  ISURF  - 1  for  surface  burst,  2  for  air  burst 
c  ZLOG  -  logarithm  (base  10)  of  scaled  range 
c 

parameter  (ns=9,  nf=7) 

real  csurf(0:ns),  cfree(0:nf) 

data  csurf /-0.173607601251,  1.35706496258, 

&  0.052492798645,  -0.196563954086, 

&  -0.0601770052288,  0.0696360270891, 

&  0.0215297490092,  -0.0161658930785, 

&  -0.00232531970294, 0.00147752067524/ 

data  cfree/-0.0423733936826, 1.36456871214, 
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& 

& 

& 


-0.0570035692784,-0. 1 82832224796, 
0.01 18851436014,  0.0432648687627, 
-0.0007997367834,-0.00436073555033/ 


c 

if(ISURF.eq.l)then 

u  =  -0.755684472698  +  1.37784223635*ZLOG 
conwep_tarr  =  csurf(ns) 
do  10  i  =  ns-1,0,-1 

10  conwep_tarr  =  conwep_tarr*u  +  csurf(i) 
else 

u  =  -0.80501734056  +  1.37407043777*ZLOG 
conwep_tarr  =  cfree(nf) 
do  20  i  =  nf-l,0,-l 

20  conwep_tarr  =  conwep_tarr*u  +  cfree(i) 
endif 

conwep_tarr  =  10.**conwep_tarr 

return 

end 


real  function  conwep_tdur(ISURF,ZLOG) 
c 

c  purpose:  find  the  scaled  duration  for  the  detonation  of  a  1  lb 
c  equivalent  hit  charge,  equations  are  from  BRL  technical 
c  report  ARBRL-TR-02555.  duration  is  returned  in 

c  msec/lb**(l/3). 

c 

c  description  of  variables: 
c  ISURF  - 1  for  surface  burst,  2  for  air  burst 
c  ZLOG  -  logarithm  (base  10)  of  scaled  range 
c 

parameter  (ns  1  =5  ,ns2=8  ,ns3=5  ,nf  1  =8  ,nf2=8  ,nf3=7) 
real  csurfl(0:nsl),  csurf2(0:ns2),  csurf3(0:ns3) 
real  cffeel(0:nfl),  cfree2(0:nf2),  cffee3(0:nf3) 
data  csurfl  /  -0.728671776005,  0.130143717675, 

&  0.134872511954,  0.0391574276906, 

&  -0.00475933664702,-0.00428144598008  / 

datacsurf2/  0.20096507334,  -0.0297944268976, 

&  0.030632954288,  0.0183405574086, 

&  -0.017396466621 1,  -0.00106321963633, 

&  0.00562060030977, 0.0001618217499, 

&  -0.0006860188944/ 

datacsurf3  /  0.572462469964,  0.0933035304009, 
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&  -0.0005849420883,  -0.00226884995013, 

&  -0.00295908591505, 0.00148029868929  / 

data  cfreel  /  -0.801052722864,  0.164953518069, 
&  0.127788499497,  0.00291430135946, 

&  0.00187957449227, 0.0173413962543, 

&  0.00269739758043,-0.00361976502798, 

&  -0.00100926577934/ 

datacfree2/  0.115874238335,  -0.0297944268969, 
&  0.0306329542941,  0.018340557407, 

&  -0.0173964666286,  -0.00106321963576, 

&  0.0056206003128,  0.0001618217499, 

&  -0.0006860188944/ 

datacfree3/  0.50659210403,  0.0967031995552, 
&  -0.00801302059667, 0.00482705779732, 

&  0.00187587272287,-0.00246738509321, 

&  -0.000841116668,  0.0006193291052  / 

c 

c  surface  burst 

if  (ISURF.eq.l)  then 
if(ZLOG  .It.  -.34 )  then 
conwep_tdur  =  -.725 
elseif(ZLOG  .It.  0.4048337)  then 
u  =  -0.1790217052  +  5.25099193925*ZLOG 
conwep_tdur  =  csuifl(nsl) 
do  11  i  =  nsl-1,0,-1 

1 1  conwep_tdur  =  conwep_tdur*u  +  csurfl(i) 
elseif(ZLOG  .It.  0.845098)  then 

u  =  -5.85909812338  +  9.299628861  l*ZLOG 
conwep_tdur  =  csurf2(ns2) 
do  12  i  =  ns2-l,0,-l 

12  conwep_tdur  =  conwep_tdur*u  +  csurf2(i) 
else 

u  =  -4.92699491 141  +  3.4634974557  l*ZLOG 
conwep_tdur  =  csurf3(ns3) 
do  13  i  =  ns3-l,0,-l 

13  conwep_tdur  =  conwep_tdur*u  +  csuif3(i) 
endif 

else 

c  air  burst 

if(ZLOG  .It.  -.34 )  then 
conwep_tdur  =  -.824 
elseif(ZLOG  .It.  0.350248)  then 
u  =  0.209440059933  +  5.11588554305*ZLOG 
conwep_tdur  =  cffeel(nfl) 
do21  i  =  nfl-l,0,-l 
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2 1  conwep_tdur  =  conwep_tdur*u  +  cfree  1  (i) 
elseif(ZLOG  .It.  0.7596678)  then 

u  =  -5.06778493835  +  9.29962886  ll*ZLOG 
conwep_tdur  =  cfree2(nf2) 
do  22  i  =  nf2-l,0,-l 

22  conwep_tdur  =  conwep_tdur*u  +  cfree2(i) 
else 

u  =  -4.39590184126  +  3.1524725264*ZLOG 
conwep_tdur  =  cfree3(nf3) 

do  23i  =  nf3-l,0,-l 

23  conwep_tdur  =  conwep_tdur*u  +  cfree3(i) 
endif 

endif 

conwep_tdur=  10.**conwep_tdur 

return 

end 


real  function  conwep_ximpr(isurf,zlog) 
c 

c  purpose:  find  the  scaled  normally  reflected  impulse  from  the 
c  detonation  of  a  1  lb  equivalent  mt  charge,  equations 
c  are  from  BRL  technical  report  ARBRL-TR-02555.  impulse 
c  is  returned  in  units  of  psi*msec/lb**(l/3) 

c 

c  description  of  variables: 
c  isurf  - 1  for  surface  burst,  2  for  air  burst 
c  zlog  -  logarithm  (base  10)  of  scaled  range 
c 

parameter  (ns=3,  nf=3) 

real  csurf(0:ns),  cfree(0:nf) 

datacsurf/  1.75291677799,  -0.949516092853, 

&  0.1 121361 18689,  -0.0250659183287  / 

data  cfree  /  1.60579280091,  -0.903118886091, 

&  0.101771877942,  -0.0242139751146/ 

c 

if  (isurf.eq.l)  then 

u  =  -0.781951689212  +  1.33422049854*zlog 
ximpr  =  csurf(ns) 
do  10  i  =  ns-1,0,-1 
10  ximpr  =  ximpr*u  +  csurf(i) 
else 
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u  =  -0.757659920369  +  1. 3788299601 8*zlog 
ximpr  =  cfree(nf) 
do  20  i  =  nf-1,0,-1 
20  ximpr  =  ximpr*u  -f-  cfree(i) 
endif 

conwep_ximpr  =  10.**ximpr 

return 

end 


real  function  conwep_ximps(isurf,zlog) 
c 

c  purpose:  find  the  scaled  incident  impulse  for  the  detonation  of 
c  a  1  lb  equivalent  mt  charge,  equations  are  from  brl 
c  technical  report  arbrl-tr-02555.  impulse  is  returned  in 
c  units  of  psi*msec/lb* *(1/3). 
c 

c  description  of  variables: 
c  isurf  - 1  for  surface  burst,  2  for  air  burst 
c  zlog  -  logarithm  (base  10)  of  scaled  range 
c 

parameter  (ns  1  =4,ns2=7,nf  1  =4,nf2=8) 

real  csurf  1  (0:ns  1 )  ,csurf2(0:ns2),cfree  1  (0:nf  1  ),cffee2(0:nf2) 

data  csurfl  /  1.57159240621,  -0.502992763686, 

&  0.171335645235, 0.0450176963051, 

&  -0.0118964626402/ 

data  csurfZ  /  0.719852655584,  -0.384519026965, 

&  -0.0280816706301,  0.00595798753822, 

&  0.014544526107,  -0.00663289334734, 

&  -0.00284189327204,  0.0013644816227  / 

data  cfreel  /  1.43534136453,  -0.443749377691, 

&  0.168825414684,0.0348138030308, 

&  -0.010435192824/ 

data  cffee2  /  0.599008468099,  -0.40463292088, 

&  -0.0142721946082,  0.00912366316617, 

&  -0.0006750681404,  -0.00800863718901, 

&  0.00314819515931, 0.00152044783382, 

&  -0.0007470265899 / 

c 

if  (isurf.eq.l)  then 
if(zlog  .It.  0.382017)  then 
u  =  0.832468843425  -h  3.0760329666*zlog 
conwep_ximps  =  csurfl  (ns  1) 
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do  11  i  =  nsl-1,0,-1 

1 1  conwep_ximps  =  conwep_ximps*u  +  csuifl(i) 
else 

u  =  -2.91358616806  +  2.40697745406*zlog 
conwep_ximps  =  csurf2(ns2) 
do  12  i  =  ns2-l,0,-l 

12  conwep_ximps  =  conwep_ximps*u  +  csurf2(i) 
endif 

else 

if(zlog  .It.  0.30103)  then 
u  =  1.04504877747  3.24299066475*zlog 
conwep_ximps  =  cfreel(nfl) 
do  21  i  =  nfl-1,0,-1 

21  conwep_ximps  =  conwep_ximps*u  -i-  cfreel(i) 
else 

u  =  -2.67912519532  2.3062923 1803*zlog 
conwep_ximps  =  cfree2(nf2) 
do  22i  =  nf2-l,0,-l 

22  conwep_ximps  =  conwep_ximps*u  -i-  cfree2(i) 
endif 

endif 

conwep_ximps  =  10.**conwep_ximps 

return 

end 
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APPENDIX  B: 

MODMCATIONS  TO  THE  DYNA2D  CODE 
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Intentionally  left  blank. 
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1.  In  subroutine  dynai,  change 


IF  (IBRODE.NE.O)  CALL  BRODEI 
to 

IF  (IBRODE.EQ.1)  CALL  BRODEI 
IF  aBRODE.EQ.2)  CALL  conwep_BLASTI 


2.  In  subroutine  fe2a,  change 
IF(LCC.GT.O)  GO  TO  15 
IF(LCC.EQ.O)  THEN 
XCTR=.5*(R(IR)+R(IL)) 

YCTR=0. 

ZCTR=.5*(Z(IR)+Z(IL)) 

CALL  SBRODE(XCTR,YCTR,ZCTR,TT,P,NPC,F) 
GO  TO  30 
ENDIF 
ISHEAR=1 
LCC-LCC 

15  IF(T(N).NE.0.0)  GO  TO  20 
to 

IF(LCC.GT.O)  GO  TO  15 
if(ibrode.eq.O)then 
ISHEAR=1 
LCC=-LCC 
else 

IF(LCC.LE.O)  THEN 
IF(LCC.EQ.-1)  THEN 
XCTR=.5*(R(IR)+R(IL)) 

YCTR=0. 

ZCTR=.5*(Z(IR)+Z(IL)) 

CALL  SBRODE(XCTR,YCTR,ZCTR,TT,P,NPC,F) 
ELSEIF  (LCC.EQ.-2)  THEN 
width=sqrt((r(ir)-r(il))**2  +  (z(ir)-z(il))**2) 
call  conwep_blast(  R(IR),  R(IL),  R(IL),  R(IR), 


& 

-width,  -width,  width,  width. 

& 

Z(IR),  Z(IL),  Z(IL),  Z(IR), 

& 

PMULT(1,N),  TT,F) 

ENDIF 

endif 

GO  TO  30 

IF(T(N).NE.0.0)  GO  TO  20 
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Intentionally  left  blank. 
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APPENDIX  C: 

MODMCATIONS  TO  THE  DYNA3D  CODE 
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Intentionally  left  blank. 
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9?  9? 


1.  In  subroutine  dynai,  change 


IF  (IBRODE.NE.O)  CALL  BRODEI 
to 

IF  (IBRODE.EQ.l)  CALL  BRODEI 
IF  (IBRODE.EQ.2)  CALL  conwep_BLASTI 

2.  In  subroutine  load,  change 

IF  (LC(l,IP+I).LT.O)  THEN 

CALLSBRODE(XCTR(I),YCTR(I),ZCTR(I),TT,P,NPC,F(I)) 

ENDIF 

to 

IF  (LC(1,IP+I).EQ.-1)  THEN 

CALLSBRODE(XCTR(I),YCTR(I),ZCTR(I),TT,P,NPC,F(I)) 
ELSE  IF  (LC(l,IP+I).EQ.-2)  THEN 
call  conwep_blast(  XXI  1(I),XX12(I),XX13(I),XX14(I), 
XX21(I),XX22(I),XX23(I),XX24(I), 

XX3 1  (I),XX32(I),XX33(I),XX34(I), 
pmult(l,ip+i),TT,F(I)) 


ENDIF 
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Intentionally  left  blank. 


40 


APPENDIX  D; 

MODMCATIONS  TO  THE  MAZE  CODE 


41 


Intentionally  left  blank. 
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1.  Insert  the  following  in  block  data: 


COMMON /BRODCM/  biodcmv(33),  ibrode 
DATAibrode/0/ 

2.  Insert  the  following  at  the  appropriate  places  in  subroutine 
meshol2: 


COMMON  /BRODCM/  brodcmv(33),  ibrode 

IF  (ICMl.EQ.'BROD'.or.icmLeq.’BLAS')  GO  TO  815 

815  call  brodin 
go  to  20 


3.  Add  subroutine  brodin: 


SUBROUTINE  BRODIN 

COMMON  /BRODCM/  YLDTI,XBO,YBO,ZBO,TBO,CL,CT,CP,IYOPT,RANGE(5), 
.COEFF(8),GFUNC(7),IUNIT,ISURF, CM, ibrode 
CHARACTER*4  ICMD 
CHARACTER*!  IPMT 
DATA  IPMT/'.y 
IBRODE=l 
10  CONTINUE 

CALL  GETSYM(1,IPMT,ICMD,4) 

IF(ICMD.NE.MODE')GO  TO  15 
CALL  GETNUM(1,IPMT,CRUD, IBRODE) 

GO  TO  900 
15  CONTINUE 

IF(ICMD.NE.  YLD'.and.icmd.ne.'WEIG')GO  TO  20 
CALL  GETNUM(1,IPMT,YLD,ICRUD) 

GO  TO  900 
20  CONTINUE 

IFCICMD.NE.HEIOGO  TO  30 
CALL  GETNUM(1,IPMTJI,ICRUD) 

GO  TO  900 
30  CONTINUE 

IFaCMD.EQ.’XBO')ICMD='XB0' 

IFaCMD.EQ.'YBO')ICMD=TB0’ 

IFaCMD.EQ.’ZBO’)ICMD='ZB0' 

IF(ICMD.EQ.TBO’)ICMD='TB0' 

IF(ICMD.NE.'XB0')GO  TO  40 
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CALL  GETNUM(1,IPMT,XB0,ICRUD) 

GO  TO  900 
40  CONTINUE 
IF(ICMD.NE.'YBO')GO  TO  45 
CALL  GETNUM(1,IPMT,YB0,ICRUD) 

GO  TO  900 
45  CONTINUE 
IF(ICMD.NE.'ZB0')GO  TO  50 
CALL  GETNUM(1  ,IPMT,ZB0,ICRUD) 

GO  TO  900 
50  CONTINUE 
IF(ICMD.NE.'TB0')GO  TO  60 
CALL  GETNUM(1  ,IPMT,TB0,ICRUD) 

GO  TO  900 
60  CONTINUE 
IF(ICMD.NE.'CM')GO  TO  65 
CALL  GETNUM(1,IPMT,CM,ICRUD) 

GO  TO  900 

65  IF(ICMD.NE.'CL')GO  TO  70 
CALL  GETNUM(1,IPMT,CL,ICRUD) 

GO  TO  900 
70  CONTINUE 
IF(ICMD.NE.'CT')GO  TO  80 
CALL  GETNUM(1,IPMT,CT,ICRUD) 

GO  TO  900 
80  CONTINUE 
IF(ICMD.NE.'CF)GO  TO  90 
CALL  GETNUM(1,IPMT,CP,ICRUD) 

GO  TO  900 
90  CONTINUE 

IF(ICMD.NE.'RANG')GO  TO  110 
DO  100 1=1,5 

CALL  GETNUM(1,IPMT,RANGE(I),ICRUD) 
100  CONTINUE 
IYOPT=l 
GO  TO  900 
110  CONTINUE 

IF(ICMD.NE.'COEP)GO  TO  130 
DO  120 1=1,8 

CALL  GETNUM(l,IPMT,COEFF(I),ICRUD) 
120  CONTINUE 
IYOPT=l 
GO  TO  900 
130  CONTINUE 

IF(ICMD.NE.'GFUN')GO  TO  150 
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DO  1401=1,7 

CALL  GETNUM(1,IPMT,GFUNC(I),ICRUD) 

140  CONTINUE 
IYOPT=l 
GO  TO  900 
150  CONTINUE 
IF(ICMD.NE.'IUNr)GOTO  152 
CALL  GETNUM(1,IPMT,CRUD,IUNIT) 

GO  TO  900 
152  CONTINUE 
IF(ICMD.NE.'ISUR')GO  TO  154 
CALL  GETNUM(1,IPMT,CRUD,ISURF) 

GO  TO  900 

158  IF(ICMD.NE.';')GO  TO  160 
RETURN 
160  CONTINUE 
WRITE(*,2010) 

900  CONTINUE 
GO  TO  10 

2010  FORMATC  ***  ERROR  ***  -  ILLEGAL  COMMAND  IN  ERODE  FUNCTION  INPUT') 
END 

4.  In  subroutine  mesholS,  after 

CALL  WRTIVL  (A(N4+N2-N1),A(N4),NUMNP0) 
add 

call  ndbrod 

5.  Add  subroutine  ndbrod: 

SUBROUTINE  NDBROD 

COMMON  /BRODCM/  YLD,H,XBO,YBO,ZBO,TBO,CL,CT,CP,IYOPT,RANGE(5), 
.COEFF(8),GFUNC(7),IUNIT,ISURF,cm,ibrode 

if(ibrode.eq.  1  )then 
write(12,'("$  ERODE  DATA")') 

WRITE(12,2010)YLD,H,XB0,YB0,ZB0,TB0 
WRITE(12,2020)CL,CT,CP,IYOPT 
IF(IYOPT.EQ.0)RETURN 
WRITE(  1 2,2030)RANGE 
WRITE(  1 2,2030)COEFF 
WRITE(  1 2,2030)GFUNC 

elseif(ibrode.eq.2)then 
write(12,'("$  CONWEP  BLAST  DATA")') 

45 


WRITE(12,2040)YLD,XB0,YB0,ZB0,TB0,IUNrr,ISURF 
if(iunit.eq.5)then 
WRITE(1 2,2040)cm,cl,ct,cp 
endif 


endif 

RETURN 

2010  FORMAT(6E10.3) 
2020  F0RMAT(3E1 0.3,15) 
2030  FORMAT(8E10.3) 
2040  FORMAT(5E10.3,2I5) 
END 
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APPENDIX  E: 

MODMCATIONS  TO  THE  INGRID  CODE 


47 


Intentionally  left  blank. 
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1 .  In  subroutine  dnopts,  change 

IF(IDATA.NE.'BROD')GO  TO  10 
to 

IF(IDATA.NE.'BROD'  .and.  idata.ne.'BLAS')GO  TO  10 


2.  Replace  subroutine  brodin  with 


SUBROUTINE  BRODIN 
character  *4  ianal,ibwm,iteo,nsnid 
common  /junkOO/  ianal,ibwm,iteo,nsmd 

COMMON  /JUNKOl/  GX,GY,GZ,NSTEP,DELT,TERM,DCTOL,ECTOL, 

.  SfflFT,MAXMEM,NMODES,IPRT,IPLT,ANIPl  ,ANIP2,NLOADC,RX,RY,RZ, 
.NBSR,NBEI,NIBSR,MSRF,ISBRF,PRTI,PLTI,IBRODE,NDETP, 
.TSSF,NM0MDP,NVELBC,AITSS,IRIGID,NUMNPB,NUMEPB,NUM1PB,NUM2PB, 
.NUM3PB,NUM4PB,ISSTRN,IHYDRO,ISHLNRM,ISHLTCK,ISHLFRM 
COMMON  /BRODCM/  YLD,H,XBO,YBO,ZBO,TBO,CL,CT,CP,IYOPT,RANGE(5), 

.  COEFF(8),GFUNC(7),IUNIT,ISURF,CM,IJUNK 
CHARACTER*4  ICMD 
CHARACTER*!  IPMT 
DATA  IPMT/'.V 
IBRODE=l 
10  CONTINUE 

CALL  GETSYM(1, IPMT, ICMD) 

CALL  SETOPT(ICMD) 

IF(ICMD.NE.'MODE')GO  TO  15 
CALL  GETNUM(l,IPMT,CRUD,IBRODE) 

GO  TO  900 
15  CONTINUE 

IF(ICMD.NE.'YLD'.and.icmd.ne.'WEIG')GO  TO  20 
CALL  GETNUM(1,IPMT,YLD,ICRUD) 

GO  TO  900 
20  CONTINUE 

IF(ICMD.NE.'HEIG')GO  TO  30 
CALL  GETNUM(1,IPMT,H,ICRUD) 

GO  TO  900 
30  CONTINUE 

IFaCMD.EQ.'XBO')ICMD='XB0' 

IF(ICMD.EQ.'YBO')ICMD='YB0' 

IF(ICMD.EQ.'ZBO')ICMD='ZB0' 

IF(ICMD.EQ.'TBO')ICMD='TB0' 

IF(ICMD.NE.'XB0')GO  TO  40 
CALL  GETNUM(1,IPMT,XB0,ICRUD) 

GO  TO  900 
40  CONTINUE 
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IF(ICMD.NE.'YBO')GO  TO  45 
CALL  GETNUM(1,IPMT,YB0,ICRUD) 

GO  TO  900 
45  CONTINUE 
IF(ICMD.NE.'ZBO')GO  TO  50 
CALL  GETNUM(1,IPMT,ZB0,ICRUD) 

GO  TO  900 
50  CONTINUE 
IF(ICMD.NE.'TB0')GO  TO  60 
CALL  GETNUM(1  ,IPMT,TB0,ICRUD) 

GO  TO  900 
60  CONTINUE 
IF(ICMD.NE.'CM')GO  TO  65 
CALL  GETNUM(1,IPMT,CM,ICRUD) 

GO  TO  900 

65  IF(ICMD.NE.'CL')GO  TO  70 
CALL  GETNUM(1,IPMT,CL,ICRUD) 

GO  TO  900 
70  CONTINUE 
IF(ICMD.NE.'CT')GO  TO  80 
CALL  GETNUM(1,IPMT,CT,ICRUD) 

GO  TO  900 
80  CONTINUE 
IF(ICMD.NE.'CP')GO  TO  90 
CALL  GETNUM(1,IPMT,CP,ICRUD) 

GO  TO  900 
90  CONTINUE 

IF(ICMD.NE.'RANG')GO  TO  1 10 
DO  100 1=1,5 

CALL  GETNUM(1,IPMT,RANGE(I),ICRUD) 
100  CONTINUE 
IYOPT=l 
GO  TO  900 
1 10  CONTINUE 

IF(ICMD.NE.'COEF')GO  TO  130 
DO  120 1=1,8 

CALL  GETNUM(l,IPMT,COEFF(I),ICRUD) 
120  CONTINUE 
IYOPT=l 
GO  TO  900 
130  CONTINUE 

IF(ICMD.NE.'GFUN')GO  TO  150 
DO  1401=1,7 

CALL  GETNUM(1,IPMT,GFUNC(I),ICRUD) 
140  CONTINUE 
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IYOPT=l 
GO  TO  900 
150  CONTINUE 
IF(ICMD.NE.'IUNr)GO  TO  152 
CALL  GETNUM(1,IPMT,CRUD,IUNIT) 

GO  TO  900 
152  CONTINUE 
rF(ICMD.NE.'ISUR')GO  TO  154 
CALL  GETNUM(1,IPMT,CRUD,ISURF) 

GO  TO  900 
GO  TO  900 

158  IF(ICMD.NE.';')GO  TO  160 
CALL  POPMG 
RETURN 
160  CONTINUE 
go  to  900 
180  continue 
WRITE(*,2010) 

WRITE(3,2010) 

CALL  TRACE(l) 

900  CONTINUE 
CALL  POPMG 
GO  TO  10 

2010  FORMATC  ***  ERROR  ***  -  ILLEGAL  COMMAND  IN  ERODE  FUNCTION  INPUT') 
END 


3.  Replace  subroutine  ndbrod  with 


SUBROUTINE  NDBROD(IBRODE) 

COMMON  /BRODCM/  YLD,H,XB0,YB0,ZB0,TB0,CL,CT,CP,IYOPT,RANGE(5), 
.  COEFF(8),GFUNC(7),IUNIT,ISURF,cm 


if(ibrode.eq.  l)then 
write(13,'("$  ERODE  DATA")') 
WRITE(13,2010)YLD,H,XB0,YB0,ZB0,TB0 
WRITE(  1 3,2020)CL,CT,CP,IYOPT 
IF(IYOPT.EQ.0)RETURN 
WRITE(13,2030)RANGE 
WRITE(  1 3,2030)COEFF 
WRITE(1 3,2030)GFUNC 


elseif(ibrode.eq.2)then 
write(13,'("$  BLAST  DATA")') 

WRITE(  1 3,2040)  YLD,XB0,YB0,ZB0,TB0,IUNIT,ISURF 
if(iunit.eq.5)then 
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WRITE(1 3,2040)cm,cl,ct,cp 
endif 


endif 


RETURN 

2010  FORMAT(6E10.3) 
2020  F0RMAT(3E1 0.3,15) 
2030  FORMAT(8E10.3) 
2040  FORMAT(5E10.3,3I5) 
END 
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Intentionally  left  blank. 
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The  following  changes  should  be  made  to  the  DYNA2D  manual,  April  1992, 

when  the  ConWep  blast  model  is  installed. 

Page  41,  Card  7,  field  6-10  (IBRODE) 
add 

EQ,2:  ConWep  blast  model  parameters  are  defined. 

Page  157,  Pressure  and  Shear  Loads 
change 

LT.O:  Erode  function  is  used  to  determine  pressure  time  history 
to 

-1:  Erode  function  is  used  to  determine  pressure  time  history 

-2:  ConWep  blast  model  is  used  to  determine  pressure  time  history 

Page  178,  Erode  Functions 
change 

if  columns  21-25  are  blank 
to 

if  columns  21-25  (lERODE)  are  blank 
change 

Otherwise,  enter  two  cards 
to 

Otherwise,  if  IBRODE  EQ  1  enter  two  cards 
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add  pages  178.1  and  178.2: 


DYNA2D  User  Manual  ConWep  Blast  Model  Input 


If  IBRODE  EQ  2,  enter  one  or  two  cards  for  the  pertinent  ConWep 
blast  function  data. 

■I - f 

I  Card  1  I 
+ - + 

Columns  Quantity  Format 


1-10 

Equivalent  weight  of  TNT 

ElO.O 

11-20 

XO  (x-coordinate  of  blast  origin) 

ElO.O 

21-30 

YO  (x-coordinate  of  blast  origin) 

ElO.O 

31-40 

ZO  (x-coordinate  of  blast  origin) 

ElO.O 

51-60 

TO  (time  of  explosion) 

ElO.O 

61-65 

lUNTT  (problem  units) 

15 

EQ.O:  iunit=2 

EQ.l:  feet,  pounds,  seconds,  psi 

EQ.2  :  meters,  kilograms,  seconds,  Pascals 

EQ.3:  inch,  dozens  of  slugs,  seconds,  psi 

EQ.4:  centimeters,  grams,  microseconds.  Megabars 

EQ.5:  user  conversions  will  be  read  from  Card  2 

66-70  ISURF  (surface  burst  or  air  burst  flag)  15 


EQ.O:  isurf=2 
EQ.l:  surface  burst 
EQ.2:  air  burst 


If  lUNTT  is  not  5,  skip  the  rest  of  this  section.  Otherwise, 
include  Card  2: 

+ - + 

I  Card  2  I 
+ - + 

Columns  Quantity  Format 

I- 10  ULB  pounds  per  problem  length  unit  ElO.O 

II- 20  UFT  feet  per  problem  length  unit  ElO.O 

21-30  UMS  milliseconds  per  problem  length  unit  ElO.O 
31-40  UPSI  psi  per  problem  pressure  unit  ElO.O 


The  ConWep  functions  are  internally  written  in  terms  of  feet, 
milliseconds,  pounds,  and  psi.  lUNTT  supplies  the  conversion 
factors  to  DYNA2D  units. 

The  development  and  limitations  underlying  this  option  are  given 
in  (Randers-Pehrson  and  Bannister,  1997), 

(Kingery  and  Bulmash,  1984)  and  (Hyde,  1993) 
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Page  241 ,  References 
add 

30.1  Hyde,  David  W.,  "User's  Guide  for  Microcomputer  Program  ConWep, 
Applications  of  TM  5-855-1,  'Fundamentals  of  Protective  Design  for 
Conventional  Weapons'",  U.S.  Army  Corps  of  Engineers  Waterways  Experiment 
Station  Instmction  Report  SL-88-1,  Vicksburg,  MS,  April  1988,  revised 
Febmary  1993. 

and 

38.1  Kingery,  Charles  N.,  and  Gerald  Bulmash,  "Airblast  Parameters  from  TNT 
Spherical  Air  Burst  and  Hemispherical  Surface  Burst,"  U.S.  Army  Ballistic 
Research  Laboratory  Technical  Report  ARBRL-TR-02555,  Aberdeen  Proving 
Ground,  MD,  April  1984. 

Page  242,  References 

add 

47.1  Randers-Pehrson,  Glenn,  and  Kenneth  A.  Bannister,  "Airblast  Loading 
Model  for  DYNA2D  and  DYNA3D,”  U.S.  Army  Research  Laboratory  Technical 
Report  ARL-TR-1310,  Aberdeen  Proving  Ground,  MD,  March  1997. 
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APPENIX  G: 

CHANGES  TO  THE  DYNA3D  MANUAL 


59 


Intentionally  left  blank. 
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The  following  changes  should  be  made  to  the  DYNA3D  manual,  November  1993, 

when  the  ConWep  blast  model  is  installed. 

Page  44,  Card  4,  field  21-25  (IBRODE) 
add 

EQ.2:  ConWep  blast  model  parameters  are  defined. 

Page  249,  Pressure  Loads 
change 

LT.O:  Erode  function  is  used  to  determine  pressure 
to 

EQ.-l:  Erode  function  is  used  to  determine  pressure 

EQ.-2:  ConWep  blast  model  is  used  to  determine  pressure  time  history 

Page  241,  Erode  Functions 
change 

if  columns  21-25  are  blank 
to 

if  columns  21-25  (lERODE)  are  blank 
change 

Otherwise,  enter  two  cards 
to 

Otherwise,  if  lERODE  EQ  1  enter  two  cards 
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add  pages  241.1  and  241.2: 


DYNA3D  User  Manual  ConWep  Blast  Model  Input 


If  IBRODE  EQ  2,  enter  one  or  two  cards  for  the  pertinent  ConWep 
blast  function  data. 


+ - + 

I  Card  1  I 
+ - + 

Columns  Quantity  Format 


1-10 

Equivalent  weight  of  TNT 

ElO.O 

11-20 

XO  (x-coordinate  of  blast  origin) 

ElO.O 

21-30 

YO  (x-coordinate  of  blast  origin) 

ElO.O 

31-40 

ZO  (x-coordinate  of  blast  origin) 

ElO.O 

51-60 

TO  (time  of  explosion) 

ElO.O 

61-65 

lUNTT  (problem  units) 

15 

EQ.O:  iunit=2 

EQ.l:  feet,  pounds,  seconds,  psi 
EQ.2:  meters,  kilograms,  seconds,  Pascals 
EQ.3:  inch,  dozens  of  slugs,  seconds,  psi 
EQ.4:  centimeters,  grams,  microseconds.  Megabars 
EQ.5:  user  conversions  will  be  read  from  Card  2 

66-70  ISURF  (surface  burst  or  air  burst  flag)  15 

EQ.O:  isurf=2 
EQ.l:  surface  burst 
EQ.2:  air  burst 
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If  lUNTT  is  not  5,  skip  the  rest  of  this  section.  Otherwise, 
include  Card  2: 


+ - + 

I  Card  2  I 
+ - + 

Columns  Quantity  Format 


I- 10  ULB  pounds  per  problem  length  unit  ElO.O 

II- 20  UFT  feet  per  problem  length  unit  ElO.O 

21-30  UMS  milliseconds  per  problem  length  unit  ElO.O 
31-40  UPSI  psi  per  problem  pressure  unit  ElO.O 

The  ConWep  functions  are  internally  written  in  terms  of  feet, 
milliseconds,  pounds,  and  psi.  lUNTT  supplies  the  conversion 
factors  to  DYNA3D  units. 

The  development  and  limitations  underlying  this  option  are  given 
in  (Randers-Pehrson  and  Bannister,  1997), 

(Kingery  and  Bulmash,  1984)  and  (Hyde,  1993) 
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Page  351,  References 
add 

44.1  Hyde,  David  W.,  "User's  Guide  for  Microcomputer  Program  ConWep, 
Applications  of  TM  5-855-1,  'Fundamentals  of  Protective  Design  for 
Conventional  Weapons'",  U.S.  Army  Corps  of  Engineers  Waterways  Experiment 
Station  Instruction  Report  SL-88-1,  Vicksburg,  MS,  April  1988,  revised 
February  1993. 


Page  352,  References 
add 


51.1  Kingery,  Charles  N.,  and  Gerald  Bulmash,  "Airblast  Parameters  from  TNT 
Spherical  Air  Burst  and  Hemispherical  Surface  Burst,"  U.S.  Army  Balhstic 
Research  Laboratory  Technical  Report  ARBRL-TR-02555,  Aberdeen  Proving 
Ground,  MD,  April  1984. 

Page  353,  References 

add 

60.1  Randers-Pehrson,  Glenn,  and  Kenneth  A.  Bannister,  "Airblast  Loading 
Model  for  DYNA2D  and  DYNA3D,”  U.S.  Army  Research  Laboratory  Technical 
Report  ARL-TR-1310,  Aberdeen  Proving  Ground,  MD,  March  1997. 


64 


APPENDIX  H: 

CHANGES  TO  THE  MAZE  MANUAL 


65 


Intentionally  left  blank. 
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The  following  changes  should  be  made  to  the  MAZE  manual,  June  1983, 
when  the  ConWep  blast  model  is  installed. 


Replace  page  36  (This  page  intentionally  left  blank)  with 
BLAST  options  ;  (DYNA2D  only)  Activate  ConWep  model 


WEIGHT  weight  Equivalent  weight  of  TNT 
XBO  X  X-coordinate  of  point  of  explosion 

YBO  y  Y-coordinate  of  point  of  explosion 

ZBO  z  Z-coordinate  of  point  of  explosion 
TBO  time  Time-zero  of  explosion 
ISURF  isurf  type  of  burst 

EQ.l:  surface  burst 
EQ.2:  air  burst  (default) 
lUNTT  iunit  Units  flag 

EQ.l:  feet,  pounds,  seconds,  psi 

EQ.2  :  meters,  kilograms,  seconds,  Pascals  (default) 

EQ.3:  inch,  dozens  of  slugs,  seconds,  psi 
EQ.4:  centimeters,  grams,  microseconds.  Megabars 
EQ.5:  user  conversions  will  be  supplied 

Conversion  factor,  pounds  per  DYNA2D  mass  unit 
Conversion  factor,  feet  per  DYNA2D  length  unit 
Conversion  factor,  milliseconds  per  DYNA2D  time  unit 
Conversion  factor,  psi  per  DYNA2D  pressure  unit 
end  of  BLAST  option  list 


CM 

ulb 

CL 

uft 

CT 

urns 

CP 

9 

upsi 

Page  27: 

PBC 

add 

If  k  is  greater  than  1,000,000,  then  k=l,000,000-k.  For  example, 
to  apply  a  ConWep  blast  load,  use  k= 1,000,002;  -2  will  be  written  in 
the  DYNA2D  input  file.  Similarly,  use  k= 1,000,001  to  apply  a  Brode 
pressure  load. 

PBCS 

add 

If  k  is  greater  than  1,000,000,  then  k=l,000,000-k  (see  PBC,  above). 
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Intentionally  left  blank. 
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APPENDIX  I: 


CHANGES  TO  THE  INGRID  MANUAL 


69 


Intentionally  left  blank. 
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The  following  changes  should  be  made  to  the  draft  INGRID  manual, 

September  1992,  when  the  ConWep  blast  model  is  installed. 

insert  on  page  page  50 

BLAST  options  ;  Activate  ConWep  model  with  the  following  optional 
parameters.  ConWep  pressure  loads  are  applied 
to  surfaces  using  load  curve  number  -2,  the  values 
for  which  are  generated  by  the  ConWep  model. 

WEIGHT  weight  Equivalent  weight  of  TNT 
XBO  X  X-coordinate  of  point  of  explosion 

YBO  y  Y-coordinate  of  point  of  explosion 

ZBO  z  Z-coordinate  of  point  of  explosion 
TBO  time  Time-zero  of  explosion 
ISURF  isurf  type  of  burst 

EQ.l:  surface  burst 
EQ.2:  air  burst  (default) 
lUNTT  iunit  Units  flag 

EQ.l:  feet,  pounds,  seconds,  psi 

EQ.2:  meters,  kilograms,  seconds,  Pascals  (default) 

EQ.3:  inch,  dozens  of  slugs,  seconds,  psi 
EQ.4:  centimeters,  grams,  microseconds.  Megabars 
EQ.5  :  user  conversions  will  be  supplied 
CM  ulb  Conversion  factor,  pounds  per  DYNA2D  mass  unit 
CL  uft  Conversion  factor,  feet  per  DYNA2D  length  unit 
CT  urns  Conversion  factor,  milliseconds  per  DYNA2D  time  unit 
CP  upsi  Conversion  factor,  psi  per  DYNA2D  pressure  unit 
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Intentionally  left  blank. 
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APPENDIX  J: 

SAMPLE  INGRID  INPUT 


73 


Intentionally  left  blank. 
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c  Surface  blast  2.32kg  TNT  y=A5.12  cm 
DN3D. 

TERM  0.001_001  PRTI  0.000_005  PLTI  0.000_005 
itss  0.000_000_005  tssf  .000_010 


blast  model  2  weight  2.32448  xbO  0  ybO  .4572  zbO  0 
tbO  -.000_100  iunits  2  isurf  1 ; 

C  IMPULSE  GAGE  MATERIAL,  SI  units 
MAT  1  1  RO  20,000.  E  0  PR  0  ENDMAT 
c  if  the  impulse  gage  is  .Im  thick,  rho=2*10**4. 
c  mass/area  =  1,000 
c  pressure  =  acceleration  *  1 ,000 
c  impulse  =  velocity  *  1,000 

mate  1 

stab  120;  1;  2;  c  impulse  gage  (10  cm  strip) 

-.1  1.1  -.10  0  -.02.02 
pr0  2  00  2  0  -2  1  0-1  0 
b  00  000  0  101000 
end 

stab  12;  12;  12  c  this  part  just  for  illustration 

-.1.1  1.1  -1.1  1.1  -.02  .02  c  2.2  m  X  2.2  m  plate 
end 
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Intentionally  left  blank. 
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APPENDIX  K: 

SAMPLE  MAZE  INPUT 


77 


Intentionally  left  blank. 
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noplot  tv  1625 
c 

Id  1  Ip  2  0  -.1  0  .15  c  axis 
Id41p200  1.00  cz=0 
Id  6  Ip  2  0  -.02  1.0  -.02  c  table  bottom 
Id  7  Ip  2  0.01  -.1  0.01  0  c  table  edge 
Id  8  Ip  2  0  -.04  1.0  -.04  c  table  bottom 
Id  9  Ip  2  0  -.06  1.0  -.06  c  table  bottom 
Id  10  Ip  2  0  -.08  1.0  -.08  c  table  bottom 
Id  1 1  Ip  2  0  -.10  1.0  -.10  c  table  bottom 

part  1  6  7  4  1  1  5  yes  c  table/impulse  gage 

part  1117  10  2  1  10  yes  c  extra  nodes  to  get  pbc/node  ratio  down  to  1/5 

c  if  the  impulse  gage  is  .02m  thick,  rho=10**6. 
c  mass/area  =  10,000 
c  pressure  =  acceleration  *  10,000 
c  impulse  =  velocity*  10,000 

assm 

title 

CONWEP  BLAST  MODEL,  HUFHNGTON  AND  EWING  TEST  T89  W=896  Z=.06322 
TERM  .000.300  PRTI  .000.0025  PLTI  .000.0025 
rcon  0 

c  put  pressure  boundary  condition  on  top  of  table 
p  1  bpbcs  3  1000002  1  1 

c  sphere  weight  896  grams  pentolite  *  .81/.7  ===  1.036  kilograms  TNT 
blast  model  2  weight  1.036  xbO  0  ybO  0 
zbO  .06322  tbO  -.000.006  iunits  2  isurf  2  ; 

c  enforce  reasonable  timestep 
itss  0.000.000.005  tssf  .000.001 

wbcd  dyna2d 

C  IMPULSE  GAGE  MATERIAL  (TABLE) 

MAT  1  1  RO  1,000,000.  E  0  PR  0  ENDMAT 

MAT  2  1  RO  1.  E  0  PR  0  ENDMAT 
t 
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Intentionally  left  blank. 
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NO.  OF 

COPffiS  ORGANTZATTON 

2  DEFNS  TECHNICAL  INFO  CTR 
ATTN  DTIC  DBA 

8725  JOHN  J  KINGMAN  RD 
STE  0944 

FT  BELVOIR  VA  22060-6218 

1  HQDA 

DAMOFDQ 

ATTN  DENNIS  SCHMIDT 
400  ARMY  PENTAGON 
WASHINGTON  DC  20310-0460 

1  US  MILITARY  ACADEMY 

MATH  SCI  CTR  OF  EXCELLENCE 
DEPT  OF  MATHEMATICAL  SCI 
ATTN  MDN  A  MAJ  DON  ENGEN 
THAYER  HALL 
WEST  POINT  NY  10996-1786 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
ATTNAMSRLCSALTP 
2800  POWDER  MILL  RD 
ADELPHI MD  20783-1 145 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
ATTN  AMSRL  CS  AL  TA 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
ATTN  AMSRL  Cl  LL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 


ABERDEEN  PROVING  GROUND 

2  DIRUSARL 

ATTN  AMSRL  Cl  LP  (305) 
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NO.  OF 

COPIES  ORGANIZATION 

6  CDR  US  ARMY  ARDEC 

ATTN  AMSTA  AR  AEE  WW 

E BAKER 

CCHIN 

RFONG 

J  PEARSON 

J  WALSH 

TECH  LIB 

PICATINNY  ARSNL  NJ  07806-5000 

1  CDR  US  ARMY  ARDEC 

ATTN  AMSTA  AR  AET  M  TECH  LIB 
PICATINNY  ARSNL  NJ  07806-5000 

1  CDR  US  ARMY  ARDEC 
ATTN  AMSTA  ARFS 
E  ANDRICOPOULOS 
PICATINNY  ARSNL  NJ  07806-5000 

3  CDR  US  ARMY  MERDEC 
ATTN  AMSME  RD  ST  WF 
L CRAFT 
D  LOVELACE 
M  SCHEXNAYDER 
REDSTONE  ARSNL  AL  35898-5250 

1  DIR  US  ARO  WASH 

ATTN  AMXRO  W  K  A  BANNISTER 
RM8N31 

5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

3  DIR  US  ARO 
ATTN  J  CHANDRA 
KIYER 

TECH  LIB 
POBOX  12211 

RSCH  TRIANGLE  PK  NC  27709-2211 

4  CDR  US  ARMY  COE 
ATTNJBALSARA 
T  BLEVINS 

P  PAPIDOS 
RNAMBURU 
3909  HALLS  FERRY  RD 
VICKSBURG  MS  39180-6199 


NO.  OF 

COPIES  ORGANIZATION 

3  DIR  US  ARMY  DARPA 
ATTN  J  RICHARDSON 
TECH  INFO 
B  WILCOX 
3701  N  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 

2  CDR  US  ARMY  TACOM 

ATTN  AMSTA  RSK  J  THOMPSON 
S  GOODMAN 
WARREN  MI  48397-5000 

2  DIRNRL 

ATTNJANEMES 
A  E  WILLIAMS  CODE  6684 
4555  OVERLOOK  AVE  SW 
WASH  DC  20375 

2  CDRNSWC 

ATTN  W  H  HOLT  CODE  G22 
WMOCK 

17320  DAHLGREN  RD 
DAHLGREN  VA  22448-5000 

8  CDR  NSWC 

ATTN  C  S  COFFEY 
RKGARRETTJR 
HMAIRR12 
J  MCKIRGAN 
BPARK 
D  G  TASKER 
TECH  LIB 
FZERILLI 

10901  NEW  HAMPSHIRE  AVE 
SILVER  SPRING  MD  20903-5000 

2  CDRNWC 

ATTN  T  GILL  CODE  3261 
TECH  LIB 

CHINA  LAKE  CA  93555-6001 

1  NAVAL  POST  GRAD  SCHL 

ATTN  J  STERNBERG  CODE  73 
MONTEREY  CA  93943 

1  USAF  PHILLIPS  LAB 

ATTN  PL  WSCD  F  ALLAHDADI 
KIRTLAND  AFB  NM  87185 


1  USAFWAL 

ATTN  T  NICHOLAS 
WRIGHT  PAT  AFB  OH  45433 
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NO.  OF 

COPIES 

4 


13 


12 


NO.  OF 


ORGANIZATION 

COPIES 

ORGANIZATION 

USAFWL 

9 

DIR  SNL 

ATTN  MNMW  W  COOK 

ATTN  R  M  BRANNON  DFV  1432 

J FOSTER 

L  C  CHHABILDAS  MS  0821 

M  NIXON 

M  FORRESTAL  DIV  1551 

TECH  LIB 

E  S  HERTEL  JR  MS  0819 

EGLIN  AFB  FL  32542-5434 

M  KIPP  DIV  1533 
JMMCGLAUNMS0819 

DIRLANL 

S  A  SILLING 

ATTN  T  F  ADAMS  F663 

T  TRUCANO  MS  0819 

J  BOLSTAD  G787 

P  YARRINGTON  DIV  1533 

J  CHAPYAK 

PO  BOX  5800 

R  DAVIDSON  K557 

EFERM 

ALBUQUERQUE  NM  87185-0307 

PFOLLANSBEEF663 

4 

lAT 

G  T  GRAY  IHB  295 

UNIV  OF  TX  AUSTIN 

KHOLIANB295 

ATTN  SJ  BLESS 

LHULL 

HD  FAIR 

J  JOHNSON  F663 

TMKEHNE 

DAMANDELLF663 

MJNORMANDLA 

P  MAUDLIN 

4030  2  W  BRAKER  LN 

TECH  LIB 

PO  BOX  1663 

AUSTIN  TX  78759-5329 

LOS  ALAMOS  NM  87454 

4 

SOUTHWEST  RSCH  INST 
DEPTOFMECHSCI 

DIRLLNL 

ATTN  C  ANDERSON 

ATTN  B  R  BOWMAN  L122 

D  LITTLEFIELD 

R  COUCH  L35 

SMULLIN 

M  FINGER  L38 

J  WALKER 

WHGOURDESr 

8500  CULEBRA  RD 

G  GOUDREAU 

PO  DRAWER  28510 

C HOOVER 

D  LASSILA  L342 

SAN  ANTONIO  TX  78284 

PRABOIN 

1 

AEROJET  ELECTRO  SYS  CO 

JEREAUGHL290 

ATTN  WARHEAD  SYSTEMS 

MJ  MURPHY 

J  CARLEONE 

TECH  LIB 

PO  BOX  296 

RE  TIPTON  L35 

PO  BOX  1663 

AZUSA  CA  91702 

LOS  ALAMOS  NM  87545 

3 

ALLIANT  TECHSYS  INC 

ATTN  S  BEISSEL 

DIRSNL 

T  HOLMQUIST  MNl  1  2720 

ATTNDBAMMANN 

RSTRYK 

MCHIESA 

G  R  JOHNSON  MNl  1  2925 

LIVERMORE  CA  94550 

7225  NORTHLAND  DR 
BROOKLYN  PK  MN  55428 

1 

ALME  AND  ASSOC 
ATTNMLALME 

102  STEVENS  FOREST  PROF 

9650  SANTIAGO  RD 

COLUMBIA  MD  21045 
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NO.  OF 

COPIES  ORGANIZATION 

1  APPLIED  RSCH  ASSOC  INC 
ATTN  J  D  YATTEAU 

5941  S  MIDDLEFIELD  RD  STE  100 
LITTLETON  CO  80123 

2  APPLIED  RSCH  ASSOC  INC 
ATTN  T  C  CARNEY 
FMAESTAS 

4300  SAN  MATEO  BLVD  SE 
STE  A220 

ALBUQUERQUE  NM  87110 

1  BRIGS  CO 
ATTNJEBACKOFEN 
2668  PETERSBOROUGH  ST 
HERNDON  VA  22071-2443 

3  DYNAEASTCORP 
ATTN  P  C  CHOU 

R  CICCARELLI 
WFLIS 

3620  HORIZON  DRIVE 
KING  OF  PRUSSIA  PA  19406 

2  GEN  RSCH  CORP 
ATTN  A  CHARTERS 
TMENNA 

5383  HOLLISTER  AVE 
SANTA  BARBARA  CA  93 1 1 1 

1  IRA  INC 

ATTN  D  ORPHAL 
4450  BLACK  AVE  STE  E 
PLEASANTON  CA  94566 

1  KAMAN  SCIENCES  CORP 
ATTN  J  S  WILBECK 
7600  BLVD  S  STE  208 
HUNTSVILLE  AL  35802 

1  KAMAN  SCIENCES  CORP 
ATTNNARI 
PO  BOX  7463 

COLORADO  SPRINGS  CO  80933-7463 

1  D  R  KENNEDY  &  ASSOC  INC 

ATTN  D  KENNEDY 
PO  BOX  4003 

MOUNTAIN  VIEW  CA  94040 


NO.  OF 

COPIES  ORGANIZATION 

1  KERLEY  PUB  SVC 
ATTNGIKERLEY 
PO  BOX  13835 

ALBUQUERQUE  NM  87192-3835 

3  LIVERMORE  SOFTWARE 
TECH  CORP 
ATTN  J  O  HALLQUIST 
B  MAKER 
D  STILLMAN 
2876  WAVERLY  WAY 
LIVERMORE  CA  94550 

1  ORLANDO  TECH  EMC 
ATTN  D  A  MATUSKA 
PO  BOX  855 
SHALLMARFL  32579 

2  SRI  INTERNATIONAL 
ATTN  D  CURRAN 

L  SEAMAN 

333  RAVENSWOOD  AVE 
MENLO  PARK  CA  94025 

1  ZERNOW  TECH  SVC  INC 

ATTNLZERNOW 
425  W  BONITA  AVE  STE  208 
SAN  DIMAS  CA  91773 


ABERDEEN  PROVING  GROUND 

1  DIR  AMSAA 

ATTN  R  THOMPSON 

49  DIRUSARL 

ATTN  AMSRL  SL  B,  P  DEITZ  (328) 

AMSRL  SL  BC,  J  T  KLOPCIC  (328) 
AMSRL  SL  BV, 

R  SAUCIER  (247) 

R  SHNIDMAN  (247) 
JRSTROBEL(247) 

AMSRL  WM,  D  ECCLESHALL 
AMSRL  WM  MA  S,  J  BEATTY 
SCHOU 
JDANDEKAR 
D  J  GROVE 
ARAJENDRAN 
T  WEERASOORIYA 
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■CQEIES 


AMSRLWMPD, 

G  GAZONAS 
D  HOPKINS 
S  WILKERSON 
AMSRL  WM  T,  W  MORRISON 
AMSRL  WM  TA, 
WBRUCHEYJR 
GFILBEYJR 
W  GILLICH 
W  GOOCH  JR 
Y HUANG 
H  MEYER  JR 
ERAPACKIJR. 

AMSRL  WMTB, 

RFREY 

J  STARKENBERG 
AMSRL  WMTC, 

TBJERKE 
R  COATES 
WSDEROSSET 
F  GRACE 
K  KIMSEY 
MLAMPSON 
L  MAGNESS 
D  SCHEFFLER 
G  SILSBY 
W  WALTERS 
AMSRL  WM  TD, 

RL  BITTING 
AM  DIETRICH  JR 
TFARRAND 
K FRANK 
N  GNIAZDOWSKI 
F  GREGORY 
P  W  KINGMAN 
MRAFTENBERG 
G  RANDERS-PEHRSON 
M  SCHEIDLER 
S  SCHOENFELD 
S  SEGLETES 
J  WALTER  JR 
T  W  WRIGHT 
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